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Abstract 

We apply a truncated set of dynamical equations of motion for connected 
equal-time Green functions up to the 4-point level to the investigation of spon¬ 
taneous ground state symmetry breaking in quantum held theory. Within 

our momentum space discretization we obtain a second order phase transition 
as soon as the connected 3-point function is included. However, an additional 
inclusion of the connected 4-point function still shows a signihcant inhuence on 
the shape of the effective potential and the critical coupling. 


1 Introduction 

In the last years the scalar quantum held theory with self-interaction has become 
an important theoretical laboratory for probing the power of non-perturbative methods 
before applying them to more involved theories as e.g. QCD. This is because $‘^-theory 
allows to study different steps of the renormalization procedure by going from (1 + 1) 
[l]-[8], over (2+1) [6]-[10] to (3+1) space-time dimensions [6, 7, 10]-[13]. 

In a previous paper [14] we investigated spontaneous ground state symmetry break¬ 
ing in the $‘^-theory in (1 + 1) dimensions by means of a closed set of equations of motion 
for connected equal-time Green functions which was derived from the coupled equa¬ 
tions of motion for the full Green functions. The necessary truncation of the originally 
inhnite set of equations was achieved by using the cluster expansions for n-point Green 
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functions, i.e. their decomposition into sums of products of connected Green functions 
(correlation functions). Within this approach we found that an inclusion of the 3-point 
function is crucial for obtaining a second order phase transition from the symmetric 
(($) = 0) to the symmetry broken (($) ^ 0) phase. This has to be the case according 
to the Simon-Grifhth theorem [15] and is in line with earlier investigations going beyond 
the Gaussian effective potential (GEP) [2]-[5]. Though the method was quite successful 
in its application to the (l + l)-dimensional theory it is not obvious what happens in 
(2+1) dimensions - which is the subject of the present paper - where we have to deal 
with a different topological structure^ leading to an additional mass counterterm in the 
renormalization procedure. 

The paper is organized as follows: In section 2 we introduce the connected Green 
function approach to the $‘^-theory in (2+1) dimensions. Section 3 contains the evalu¬ 
ation of the effective potential at zero temperature within various limits while section 
4 is devoted to a summary. The complete set of dynamical equations in conhguration 
space is presented in the appendix. 


2 Correlation dynamics for <l> 2 ^^-theory 


The Lagrangian density and the corresponding Hamiltonian are given by 


£ 

E 


— J (P X 


- 

TT^ + (V^)^ + 


•} 


( 1 ) 

( 2 ) 


respectively, with tt = and = m^ + hm^, where niQ is the bare mass, m the renor¬ 
malized mass, and 6m^ is the mass counterterm. As in the (1 + 1) dimensional case the 
$ 2 - 1-1 superrenormalizable and no coupling or wave function renormalization 

occurs. However, in addition to the divergent tadpole graph (hg.la) in (1 + 1) dimen¬ 
sions, in (2+1) dimensions we have to subtract the divergent contribution due to the 
setting sun diagram (hg.lb). The corresponding divergent self-energies are given by 


^tadpole = 3A 


(Pp 

(27r)^ 


i 

p2 _ jjpi _j_ 


(3) 


(Pp f (Pq 


j.™t(^2,^2) ^ -6l\^ J J 


(27r)3 J (27r)^ 


p2 _ jyiP _|_ 


q2 _ jyiP _|_ 


[k — p — qY — rrP + ie 


(4) 


^In the symmetry broken phase the two minima of the effective potential are unitary equivalent in 
^t+i -theory (which can be seen by the non-vanishing vacuum tunneling amplitude), while this is not 
the case in (2-1-1) dimensions (see e.g. [6, 16]). 
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a) 




Figure 1: The primitively divergent self-energy graphs in the $ 2 - 1-1 field theory; tadpole 
(a), setting sun (b). 


j^tadpoie ^^ 2 '^ is linearly divergent and has the peculiarity of being independent of the 
external momentum, while the logarithmically divergent part of m^) is given 

by the hrst term of the corresponding Taylor expansion in the external momentum 
around = 0: 


■Y^sunset / ^2 


(F,m") = ^ 

n=0 

S sunset/Ci 2\ , s^sunset / i 2 

(U, m j + L (A: , m 


2\n 


^ Y-sunsetn 2 2\ I ^ / 

dnk2^ {k,ni}\k2=0 


(5) 


Since m^) is hnite with m^) 

zero-momentum-subtraction-scheme is given by 

+ Sml^set = 


= 0 the mass counterterm in the 
(0,m^). (6) 


The subtraction of divergencies dehned in this way is completely equivalent to the 
familiar BPH or BPHZ^ renormalization scheme [17, 18]. 

Before subtracting the isolated divergencies within the method applied here, we 
have to perform the integrations over dpo in and over dpo^dqo in 

The results are: 


8m 


2 

tadpole 


-3A 
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-3AA(0), 
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A(Ti — X 2 ) = 


d^p 1 

(27r)2 2ijjp 
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^Bogoliubov-Parasiuk-Hepp-Zimmermann 
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With the divergencies specihed as above, we can now derive the closed set of equations 
of motion for equal-time Green functions up to 4-point functions. 

We start by decomposing the held ^ into a classical and a quantum part according 
to 


0 = *0 + $, ($) = 0, 


( 10 ) 


where $o is a real constant. After this separation the Hamiltonian is given by 

H = ^J d^x [n^ + (V$)2 + 2{%m^ + + (^2 ^ ^ 


( 11 ) 


with TT = dt(f> = dt^ = H. Due to the inclusion of the mass counterterms in (If) the set 
of equations of motion generated by this Hamiltonian is already properly renormalized 
since the inhnite constant zero-point contributions (still contained in the Hamiltonian) 
do not enter the Heisenberg equation. 

By means of the equal-time commutation relations 




[n(f,f), n(y, t)] = 0, [$(f,f), n(y, t)] = iS^{x - y), 


( 12 ) 


and the Heisenberg equation of motion we obtain 


dt<5 = H, 


d.n = + hm^) - + (V" - - 3A$2)$ _ 3A$o$" - A^^ (13) 

for the time evolution of $ and H. In the same way we are now able to derive equations 
of motion for equal-time operator products of $ and H of the general form 

Gn{xi,...,Xn,t) = (14) 

n 

where n is the number of held operators and all time arguments are taken to be equal 
as indicated by the single time argument t occurring in Gn (e.g. for n = 2 we obtain 
equations of motion for the operator products $(xi, f)$(x 2 G )5 n(®iG)^(® 2 G )5 ^-iid 
n(xi, f)n(x2G))- After taking the expectation values of these equations we obtain an 
inhnite hierarchy of equations of motion for equal-time Green functions 

d,Gn = cMGn) = -{[Gn,H]), (15) 

I 

where the equations of motion for n-point functions couple to (n+1)- and (n+2)-point 
functions due to the and terms in the Hamiltonian (11). We note in passing 
that for the description of the time evolution of a system in a local held theory it is 
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sufficient to consider the equal-time Green functions, as long as one considers all Green 
functions containing the helds as well as their conjugate momenta. 

In order to allow for a practical application, as e.g. the evaluation of the effective 
potential, we have to truncate the inhnite set of equations. This can be achieved using 
the cluster expansion of n-point Green functions, i.e. their decomposition into sums 
of products of connected Green functions. The explicit form of the cluster expansions 
can be derived from the representation of the theory in terms of generating functionals 
of full and connected Green functions (see e.g. [14, 18, 19]). 

Due to the separation of the held into a classical and a quantum part according 
to (10) all 1-point functions in the cluster expansions can be omitted^. Since our aim 
is a closed set of equations on the 4-point level and the equations of motion for the 
4-point functions couple to the 5- and 6-point functions, we have to specify the cluster 
expansions up to sixth order, i.e. 

( 12 ) = ( 12 )„ 

(123) = (123),, 

4 

(1234) = (1234), + (1 + ^iP2n)(12),(34)„ 

?4 = 3 

5 4 5 

(12345) = (12345), + (1 + + V2n] + E E tPmiP24(12)c(345)„ 

n=3 71=3^=71+1 

6 3 

(123456) = (123456), +(1 + EEn+(123),(456), 

71 = 4 j = l 

6 5 6 

+ (1 + E 3“ ”^2ra] + E E ”^lra”^2j)(12)c(3456), 

71 = 3 71 = 3^ = 71 + 1 

4 4 

+ (1 + E ^5n)(l + E tP2n)(12),(34),(56),. (16) 

71 = 1 71 = 3 

In the cluster expansions (16) the integer numbers between the brackets stand for 
held operators $(z) or conjugate momenta 11(1) while Vi^ is the two-body permutation 
operator. We note that any application of the Vi^ operators has to ensure the original 
order of the held operators within the connected Green functions occurring on the r.h.s. 
of (16). For a representation of the cluster expansions including the 1-point functions 
the reader is referred to ref. [19]. 

^We get the same result by using the Hamiltonian (2) instead of (11) and the complete cluster 
expansions with (d) = 4’o. 
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The cluster expansions are truncated by neglecting all connected n-point Green 
functions with n > 4, i.e. (12345)c and (123456)c in (16). Inserting the truncated 
cluster expansions into the hierarchy of equations of motion (15) up to the 4-point 
functions on the l.h.s and hence up to the 6-point functions on the r.h.s. of (15) leads 
to the desired closed set of non-linear coupled equations of motion for the connected 
Green functions, i.e. a system of correlation dynamical equations'^. 

Due to the inclusion of the mass counterterms in (11) the resulting set of equations of 
motion up to the 4-point level is properly renormalized as mentioned before. However, 
normal ordering of all operator products occurring in the equations of motion with 
respect to the mode expansions containing the renormalized mass m (perturbative 
vacuum) 



shows that additional terms only arise within the connected 2-point functions, i.e. 
($(Ti)$(T2))c = (: $(Ti)$(T2) :)c + A(Ti - T 2 ), 

(n(Ti)n(T2))c = (: n(Ti)n(T2) Oc - (v|^ - m^)A(Ti - T 2 ), 


(n(Ti)$(T2))c = (: n(Ti)$(T2) :)c - - T 2 ), 

($(Ti)n(T2))c = (: <i)(Ti)n(T2) :)c + ^<^(^1 - ^ 2 ), (18) 

where A(Ti — X 2 ) is given in (9). This can be seen by a detailed analysis taking into 
account the cluster expansions (16). 

After inserting the normal ordered 2-point functions (18) into the equations of mo¬ 
tion, all terms proportional to A(0) arising from the normal ordering of ($(Ti)$(Ti))c 
cancel out due to the counterterm (7). But contrary to the theory in (1 + 1) 

dimensions, where the tadpole is the only divergent diagram [1, 9, 14], we are left with 
the contributions from the counterterm (8) which do not cancel in the nor¬ 

mal ordered equations of motion®. The resulting equations of motion for the 2-point 

■^This lengthy but straightforward derivation can most easily be performed by means of computer 
algebra, as e.g. provided by systems as MATHEMATICA. 

®Chang and Magruder [9] have shown that even in -theory one is able to construct a Rnite 
Hamiltonian by introducing a new normal ordering prescription which, however, has to depend on the 
interaction explicitly. 
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connected Green functions read: 


|{: <l.(l)<i.(2) :), = (: n(l)<i.(2) :), + (: <i>(l)n(2) :)„ (19) 

|{: n(l)<i.(2) :), = {: n(l)n(2) :). + i(l)(: <l.(l)<i.(2) + 3A<i.S)A( 12) 

-A13(: <l.(l)<i.(l) :),({: <l.(l)<i.( 2 ) :), + A(12)) + {: $=(1 )<1.(2) :)J 

-3At„{; t(l)«(l)t(2) ;)„ (20) 

^{:n(l)n( 2 ):), = i(l)(:«(l)n( 2 ) :), 

-A13(: <I>(l)<i>(l) :).(: <I>(l)n(2) :),+ {: <i.’(l)n(2) :)J 
-3A$„(:$(l)«(l)n(2) :), 

+i(2)(:n(l)<i.(2) :), 

-A13(: n(l)<i.(2) :),{: <i.(2)<l.(2) :),+ {: n(l)<l.'’(2) :)J 

-3A$„(:n(l)«(2)$(2) :)„ (21) 

with the modihed kinetic energy 

i{{) = til) - (22) 

Although normal ordering only affects the connected 2-point functions, we have written 
out the normal ordering operation : • : in all connected Green functions in order to 
have a uniform notation. Due to their length the equations of motion for the connected 
3- and 4-point functions are not explicitly stated here and presented in the appendix. 

Finally we expand the renormalized set of correlation dynamical equations of motion 
in an arbitrary orthonormal single-particle basis {V’a}, he. 

n(F) = (23) 

O' O' 

which simplihes their numerical integration and additionally allows to discretize 

in a consistent way. In order to compactify the resulting equations we introduce the 

following abbreviations: 

(a|AiA2) = j (fx V’*(®)V’Ai(®)V’A2(®) , 

(a|AiA2A3) = j (f X tfjl{x)tfjx,ix)tfjx2ix)tfjxsix) , 


iaf3 = / d'^x i/ilix) (v^ - - 3A$2 _ V’y(a^) , 
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( 24 ) 


f/.S = -3A J;{q|A,A2;3)(; ;), 

Ai A2 

The permutation operator interchanging the indices a and /3 is denoted by Vaji- Again 
we only present the equations on the 2-point level while the 3-point and 4-point func¬ 
tions follow in a straight forward manner from the appendix. 

4(: :)c = (: :), + (: :), , (25) 

at 

:)c = (: n„n^ :)c 
at 

+ ^(faA + Ua\){- ■) c + ''^{UaX ~ ^f^^sunset^ ct\ ~ 3A$Qh„A)AA/3 

A A 

-A ^ (a|AiA2A3)(: $Ai$A 2 ^A 3 ^y :)c 

AiA2A3 

-3A$o ^ (a|AiA 2 )(: ^Ai^A^^y :)c , (26) 

Ai A2 


— (: n„n^ :)c — (1 + Vap) ^(faA + Ua\){- :)c 

dt Y 

— A(1 + "Pa^) ^ (a|AiA2A3)(: $Ai4)A2^A3n^ :)c 

AiA 2 A 3 

-3A$o(l + Vap) E («AiA 2 )(: $A 3 $A 2 n )3 :)c. (27) 

Ai A 2 

The resulting set of equations of motion for the $ 2 _|_i-theory will be denoted as $ 2 _|_iC*T* 
($2-1-1 Correlation Dynamics) furtheron. It conserves the total energy and momentum 
as can be shown within the general way outlined in [19]. 

3 Numerical solution of the ^ 2 + 1 ^^ equations 

In this section we present the application of the $ 2 - 1 - 1 ^-^ equations to the evaluation 
of the effective potential of the $ 2 +i"fheory at zero temperature, which corresponds 
to the determination of the ground state energy for a given magnetization $o- Such 
a task automatically includes an investigation of the spontaneous symmetry breaking 
from the symmetric ((^) = $0 = 0 ) to the symmetry broken ((^) = $0 7 ^ 0 ) phase, i.e. 
a spontaneous breakdown of the symmetry under the transformation ^ ^ as soon 

as the coupling exceeds a critical value. 

Since we are interested not only in the effective potential generated by the full set of 
equations but additionally in the influence of the 2-, 3- and 4-point functions we intro¬ 
duce 4 limiting cases of the $ 2 - 1 - 1 ^-^ equations, denoted by $ 2 +iC*T*( 2 ), $ 2 +iC*T*( 2 , 3), 
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$ 2+1 4) and $2+iC*-D(2, 3, 4). The numbers in parentheses indicate the various 

n-point functions taken into account, e.g. in $2+iC*-D(2, 3) the 4-point functions are 
set equal to zero while $2+iC*-D(2, 3, 4) corresponds to the full set of equations. 

At this point we again have to come back to the renormalization, since we only 
obtain meaningful results if we subtract the divergencies inherent in the introduced 
limiting cases. The divergent setting sun graph can only be generated when considering 
the 4-point function while the divergent contribution due to the tadpole is included in 
all 4 limiting cases of the $ 2 + 1 ^-^ equations. Thus the mass counterterm is simply 
given by 6m^ = (as in the case of the (1 + 1) dimensional theory) as long 

as we restrict our calculations to the limits $2+iC*T*(2) and $2+iC*T*(2, 3), hence we 
have to neglect all terms containing in the equations of motion derived in 

section 2. However, in the C'i)(2,4) and $2+iC*T*(2,3,4) approximations our 
mass counterterm consists of the divergent contributions of both, the tadpole and the 
setting sun graph. 

3.1 Application to the effective potential 

In order to solve the set of equations (25) - (27) and those for the 3- and 4-point 
functions numerically, we discretize the momentum space by choosing plane waves in a 
two-dimensional box with periodic boundary conditions as a single-particle basis, i.e. 
we consider the theory on a torus. In order to select an appropriate box size for a given 
renormalized mass we have to make sure that our results converge as a function of the 
number of single-particle states included. Due to restrictions concerning computer time 
and memory space we have to compromise between a good momentum space resolution 
and the number of single-particle states involved. For practical purposes we therefore 
use a boxlength of 20 fm and include the 29 lowest lying single-particle states for a 
renormalized mass of 10 MeV. The convergence of our numerical studies with respect 
to the number of single-particle states will be demonstrated below (section 3.2). 

As discussed in [14] there is no easy access to the stationary solutions of the cor¬ 
relation dynamical equations. Hence we are only able to describe the propagation of 
the system in time for a given initial conhguration. In order to evaluate the ground 
state equal-time Green functions within the $2+1^4) approximation we thus exploit 
the Gell-Mann and Low theorem [20] by using the method of adiabatically switching 
on the coupling as was already done in [14]. Starting with the trivial exact ground 
state conhguration for A = 0 and a given hxed value of the magnetization $0 as an 
initial condition, we are then able to obtain the ground state equal-time Green func¬ 
tions for hnite A as long as we switch on the coupling slow enough to ensure that the 
system adiabatically follows the trajectory of the ground state to a good approxima¬ 
tion. The time-dependence of the coupling is chosen to be linear with A/(4m) = l3t and 
/3 = const, but actually the functional form of the time-dependence is not essential as 
long as fi is small enough to guarantee adiabaticity. The convergence of the method. 
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Figure 2: Ground state energy density as a function of the coupling for different values 
of the dimensionless vacuum magnetization; all curves (except for the GEP approxi¬ 
mation) have been evaluated time-dependently with fi = 0.05 c/fm. 


depending on /3, was discussed in detail in [14]; it works equally well in (2-|-l) dimen¬ 
sions so that we do not present the convergence of the method again. For all results 
shown in the following we have chosen fi = 0.05 c/fm, which yields a quite good com¬ 
promise between the convergence of the method and the necessary computation time. 
Furthermore, in the subsequent discussions, we have replaced the D{2) limit of 

the correlation dynamical equations by the GEP approximation, since this represents 
the corresponding stationary limit. In order to present our results we use dimensionless 
quantities, e.g. {Hjm^) jL'^ for the energy density, Xl[4:m) for the coupling, and 
for the magnetization. 

Fig.2 shows the ground state energy density obtained as described above as a func¬ 
tion of the coupling for various values of the magnetization. We observe that in all 
approximations there exists an intersection between the ground state energy density 
curve for $o = 0 and a curve generated with a hnite value of the magnetization, i.e. 
the energy minima are obtained for non-vanishing values of the magnetization when 
the coupling exceeds a critical value. Hence all limiting cases of the ^ 2 + 1 ^'^ equa¬ 
tions predict a phase transition with increasing coupling. This result is in line with 
Magruder [9], who proves the existence of a phase transition in a properly renormal- 
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ized (2+1) dimensional $‘^-theory by extending earlier arguments of Chang [1] for the 
(1 + 1) dimensional case. The proof is based on a comparison of two Lagrangians, one 
with positive and one with negative value of so that the corresponding classical 
potentials obtain their minimum for $o = 0 and $o 7^ 0, respectively. One observes 
that a strong coupling theory dehned by the Lagrangian with negative is equivalent 
to a weak coupling limit of the Lagrangian with positive and vice versa. Therefore 
the theory described by either of the two Lagrangians has to undergo a phase tran¬ 
sition when the coupling is increased far enough. This concept of equivalent or dual 
Lagrangians was also exploited in [2] in case of the (1 + 1) dimensional $‘^-theory. Ac¬ 
cording to this concept, which does not rely on any approximative method, the phase 
transition is governed by the renormalization, i.e. by the mass counterterms introduced 
above. We will come back to this point later and proceed by recognizing that we obtain 
a hrst order phase transition in case of the GEP and $2+iC*T*(2, 4) limits, while the 
^ 2 - 1-1 3) and $2_|_iGi)(2, 3, 4) approximations lead to a phase transition of second 

order. This can seen by the fact that in case of the GEP and $2_|_iGi)(2, 4) approach 
the curve with vanishing magnetization is hrst intersected by curves with considerably 
large values of the magnetization which implies that there is a discontinuity if we follow 
the trajectory of the lowest possible energy density, i.e. the ground state magnetiza¬ 
tion as a function of A has to jump from 0 to a hnite value at the critical coupling. 
This is not the case for the $2_|_iGi)(2, 3) and $2_|_iGi)(2, 3, 4) approximations, where 
all displayed curves with increasing magnetization intersect the curve with $o = 0 in 
the correct order to yield a continuous increase of the ground state magnetization as 
a function of A. But in contrast to the (1 + 1) dimensional theory we are not able to 
determine the proper order of the phase transition by simple considerations [1], since 
the Simon-Grifhth theorem [15] is not applicable here. 

In hg.3 the different effective potentials are plotted as a function of the magneti¬ 
zation $0 for some selected values of the coupling. The data in this hgure are the 
same already used in the hg.2. In case of the GEP and $2_|_iGi)(2, 4) approximations 
for couplings around the critical value there is an intermediate maximum between two 
minima located at $o = 0 and at a hnite value of the magnetization. The interme¬ 
diate maximum leads to a sudden jump of the ground state magnetization from 0 to 
a hnite value as soon as the critical value of the coupling is reached, yielding a phase 
transition of hrst order. In contrast to this the $2_|_iGi)(2, 3) and $2_|_iGi)(2, 3, 4) ap¬ 
proximations exhibit no intermediate maximum allowing for a continuous increase of 
the magnetization. 

In addition to the order of the phase transition the different shapes of the effective 
potentials in hg.3 help to specify the regions where the 2-, 3- and 4-point functions 
reveal their main importance. Obviously, the 2-point function cannot lower the energy 
density around $o = 0, while it becomes quite important for large values of $ 0 - Due to 
the symmetry under the transformation cf) ^ —cf) the 3-point function has to vanish at 
$0 = 0, thus revealing its main importance at some intermediate value of $ 0 - Within 
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Figure 3: Effective potential as a function of the vacuum magnetization for different 
values of the dimensionless coupling A/(4m); all curves (except for the GEP approxi¬ 
mation) have been evaluated time-dependently with fi = 0.05 c/fm. 


our approximation the 4-point function is the only contribution leading to 

a reduction of the energy density around $o = 0. These qualitative considerations 
illustrate why we obtain a phase transition of second order as soon as we include the 

3- point function in our calculation. The 3-point function brings down the intermedi¬ 
ate maximum (in the energy density) in a region where neither the 2-point nor the 

4- point function can reduce the energy density substantially, yielding an effective po¬ 
tential characterized by only one minimum. A similar result concerning the relative 
importance of the different correlation functions was obtained in [14]. Therefore the 
qualitative picture that the n-point functions act dominantly at lower values of the 
magnetization than the (n-l)-point functions is conhrmed again. 

In order to present the phase transitions in a more common way the ground state 
magnetization (the order parameter of the phase transition) is plotted as a function of 
the coupling in hg.4. For this representation we again use the numerical data of hgs. 2 
and 3. In case of the GEP and $2_|_iGi)(2, 4) approximations there is a discontinuity 
at the critical coupling yielding a sudden jump of the ground state magnetization from 
0 to 0.89 and 1.1, respectively. Including the 3-point function in the $2_|_iGi)(2, 3) 
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Figure 4: Ground state magnetization as a function of the dimensionless coupling for 
the different limiting cases of the D equations. 

and $2_|_iGi)(2, 3, 4) approximations we obtain a smooth increase of the ground state 
magnetization indicating a phase transition of second order. 

Furthermore, hg.4 shows that in the $2_|_iGi)(2, 3, 4) limit the phase transition is 
of ’weak’ second order, meaning that the slope of the order parameter just above the 
critical coupling is too high in order to allow for a clear distinction from a hrst order 
transition. This observation can again easily be understood within our simple consider¬ 
ations concerning the shape of the effective potential. Since the 4-point function lowers 
the energy density around $o = 0, the intermediate maximum, responsible for the hrst 
order phase transition in the GFP and $2_|_iGi)(2, 4) limit, becomes more pronounced 
when the 4-point function is included. 

Furthermore, hg.4 shows how the inclusion of the different n-point functions changes 
the critical coupling. The extracted numerical values are given in tab.f in compari¬ 
son to values obtained in [14] for the (1 + 1) dimensional $‘^-theory. Since the present 
investigation does not correspond to the continuum limit (see below) the actual quan¬ 
titative values are not of primary interest, but a comparison of the relative values of 
the different approximations for the different dimensionalities allows to extract the in- 
huence of the divergent setting sun graph, which is only included in the $2_|_iGi)(2, 4) 
and $2_|_iGi)(2, 3, 4) approximations. Obviously the inclusion of this diagram leads to 
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(1 + 1) dimensions 

(2+1) dimensions 

approximation 

(A/4/77, ^crit 

approximation 


GEPi+i 

2.568 

GEP2+1 

1.088 

^t+iCD{2,3) 

1.629 

^l^,CD{2,3) 

0.585 

$4+iGP(2,4) 

3.81 

$4^iGP(2,4) 

0.484 

$4+iGP(2,3,4) 

2.446 

$4^iGP(2,3,4) 

0.369 


Table 1: Values of the critical coupling for different approximations and space-time 
dimensions; the values in (1 + 1) dimensions, which are taken from [14], have been 
evaluated using L = 100 fm boxlength, /3 = 0.05 c/fm, m = 10 MeV, and 15 single¬ 
particle states. 


a reduction of the critical coupling as can be seen by the fact that the $2_|_iC'i)(2, 4) 
and $2_|_iC'i)(2, 3, 4) approximations have lower values of the critical coupling than the 
D{2, 3) limit in contrast to the results obtained in the (1 + 1) dimensional theory. 

3.2 Numerical convergence 

In order to get some feeling concerning the convergence of our results we show in hg.5 
the critical coupling as a function of the momentum cutoff for different values of the 
boxsize for the GEP and $2_|_iGi)(2, 3). 

One clearly sees that the momentum discretization used for the presented results 
(L=20 fm boxlength, (kjAk)max = 3 which is equivalent to 29 single-particle states) 
is quite away from the continuum limit, where (A/4m)crU = 3.078 is the analytically 
known limit for the GEP [6] (straight line), but surely convergent with respect to the 
ultra-violet momentum cutoff. We furthermore observe that reducing the size of the 
momentum cell Ak by increasing the boxlength (keeping the number of single-particle 
states hxed) yields higher values for the critical coupling while a higher value of the 
momentum cutoff - when increasing the number of basis states and keeping the size of 
the momentum cell hxed - has the opposite effect. Thus the analysis in the present work 
does not correspond to the continuum limit of the theory, but clearly demonstrates the 
relative importance of 2-, 3- and 4-point functions for the ground state energy density 
in different regimes of the coupling as well as for the order of the phase transition. 

3.3 Comparison to related approaches 

Einally we would like to discuss the results of our present work in the context of the lit¬ 
erature known to us which also deals with the application of methods to the $‘^-theory 
in (2+1) dimensions, which go beyond the GEP approximation. Stancu and Stevenson 
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1{S} 2(13} 3(29) 4 (49) 5 (81) 6(113) 


{# of states} 

Figure 5: Critical coupling as a function of the momentum cutoff (i.e. the number of 
single-particle states) for different boxsizes; the values for the $2+iC*-D(2, 3) approxi¬ 
mation have been evaluated time-dependently with fi = 0.05c//m; the horizontal line 
in the upper plot indicates the critical coupling for the GEP approximation in the 
continuum limit. 

[7] have formulated a systematic non-perturbative expansion for the effective potential 
of the $‘^-theory based on the variational GEP method, thereby recovering, as a hrst 
order result, the GEP results. But in contrast to similar approaches as e.g. presented 
in [5, 8], the variational mass parameter is not hxed to the Gaussian value if one pro¬ 
ceeds to higher orders of the expansion, which is, according to the authors, crucial for 
the expansion to yield convergent results. Going to the second order of this expansion 
(denoted by PGEP in [7]) yields a second order phase transition to the symmetry bro¬ 
ken phase at a critical coupling of 2.775 or 2.237, for approaching the phase transition 
from below or above, respectively. In contrast to this result Cea and Tedesco [8] found 
that there is no spontaneous symmetry breaking in the (2-|-l) dimensional $‘^-theory 
by evaluating the two-loop corrections to the GEP approximation. This, however, is 
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in contradiction to our result and the proof given by Magruder [9]. The discrepancies 
arise due to the different use of the gap equation in [7] and [8], i.e. Stancu allows the 
variational mass parameter to change by going to the next order of the expansion while 
Cea keeps this value hxed. Such problems do not come up in the (1 + 1) dimensional 
theory where one only has the divergent tadpole graph which can consistently be sub¬ 
tracted in the GEP scheme yielding no ambiguities in the renormalization procedure 
when going beyond the GEP approach. In (2+1) dimensions difficulties occur when 
using the GEP as a starting point for a post Gaussian approach, since the GEP ap¬ 
proximation takes into account only parts of the diagrams included in the second order 
of the loop expansion, e.g. the cactus graphs (see e.g. [1]), while the divergent setting 
sun diagram is not considered. 


4 Summary 

In this paper we have presented the application of a connected Green function approach 
(denoted by ^‘^CD) to the evaluation of the effective potential in (2+1) dimensional 
$‘^-theory and thereby to the investigation of the ground state symmetry breaking. 

By comparing 4 limiting cases of the full D equations we were able to re¬ 

construct where the 2-, 3- and 4-point functions dominantly influence the shape of the 
effective potential, which in turn allowed to give an intuitive explanation for the order 
of the phase transition obtained. The 2-point function mainly decreases the ground 
state energy for large values of the magnetization while the 3-point function exhibits 
its main influence at intermediate values. Both are not able to reduce the ground state 
energy density for vanishing magnetization ($o —^ 0) which, however, is the region 
where the 4-point function acts dominantly. 

This picture provides a simple explanation of our results, showing a hrst order phase 
transition in the approximations which do not take into account the 3-point functions 
while the phase transition is of second order as soon as the 3-point function is included. 
The emerging picture is in line with our earlier results in [14] where the D approach 
was applied to the (1 + 1) dimensional $‘^-theory. 

In addition to the divergent tadpole graph of the $‘^-theory in (1 + 1) dimensions 
a proper renormalization in the (2+1) dimensional theory has to take into account 
the divergent setting sun graph, which is of second order in the coupling. Within our 
approximation schemes we had to care about this graph as soon as the 4-point function 
was included in our calculations, yielding a large reduction of the critical coupling. 
Therefore the $2_|_iGi)(2, 4) and $2_|_iGi)(2, 3, 4) approximations lead to lower values 
of the critical coupling than the $2_|_iGi)(2, 3) limit, in contrast to the results obtained 
in (1 + 1) dimensions. This clearly shows that the phase transition is governed by the 
counterterms taken into account. Thus the proper renormalization in each space-time 
dimensionality is essential for a correct description of the phase transition while the 
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different approximations yield similar shapes for the effective potentials. 

Since the description of the phase transition by the ^'^CD equations leads to reason¬ 
able results in (1 + 1) and (2+1) dimensions and is in line with rigorous results known to 
us (e.g. the Simon-Grifhth theorem [15]), we conclude, that correlation dynamics can 
be a powerful tool for the description of low-energy phenomena in local held theories. 
These quite encouraging results in the ’$‘^-laboratory’ thus seem to justify the applica¬ 
tion of the method to SU(N) gauge theories, aiming at a non-perturbative description 
of QCD phenomena. First steps in this direction have already been performed in [19]. 
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Appendix 

Equations of motion for the connected 3- and 4- 
point functions in configuration space 

In this appendix we present the renormalized correlation dynamical equations for the 

normal ordered connected 3- and 4-point Green functions in conhguration space in 

addition to (19)-(21) using the same notation: 

f(: timvm y = {: n(l)<l.(2)«(3) y 

+ (: $(l)n(2)$(3) :),+ (: «(l)$(2)n(:i) :)„ (28) 

4(: n(i)«(2)<i.(3) y = (: n(i)n(2)«(3) y + (: n(i)<i.(2)n(3) y 

+((!){: <1>(1)*(2)<1.(3) y 

-3A[(:$(l)«(l)y.(:$(l)«(2)$(3) :), 

+ (: <1.(1)«(1)<1.(2) y((: «(1)<1.(3) y + A(13)) 

+ ({: «(1)<1.(2) y + A(12)){: <1.(1)«(1)<1.(3) y] 

-3A«„[2({; t(l)«(2) y + A(12))((; *(1)*(3) + A(13)) 

+(: <1>(1)«(1)<1>(2)«(3) yl, (29) 
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^(: n(l)n(2)$(3) :), = (: n(l)n(2)n(3) :), + (1 + Vr,)i{l){: $(l)n(2)$(3) :), 

-3A(1 + Vr 2 )[{-. :)c(: $(l)n(2)$(3) :), 

+ (:$(l)$(l)n( 2 ) :),((: $( 1 )$( 3 ) :), + A(13)) 

+ (:$(l)n(2) :),(:$(1)$(1)$(3) :),] 

-3A$o(l + n2)[2(: $(l)n(2) :),((: $(1)$(3) :), + A(13)) 

+ (:$(l)$(l)n(2)$(3) (30) 

^(: n(l)n(2)n(3) :), = (1 + Vr 2 + Vr,)i{l){-. $(l)n(2)n(3) :), 

-3A(1 + Vr 2 + VrM- ^(1)^(1) :)c(: $(l)n(2)n(3) :), 

+ (:$(l)$(l)n(2) $(l)n(3) :), 

+ (:$(l)n(2) :),(:$(l)$(l)n(3) :),] 

-3A$o(l + V ^2 + n 3 )[ 2 (: $(l)n( 2 ) $(l)n(3) :), + (: $(l)$(l)n(2)n(3) :),] 

+ ^A$o[<^(l - 2).5(1 - 3) - .5(1 - 2)A(2 - 3) + .5(1 - 3).5(2 - 3)], (31) 


^(: $(1)$(2)$(3)$(4) :), = (: n(l)$(2)$(3)$(4) :), 

+ (: $(l)n(2)$(3)$(4) :), + (: $(l)$(2)n(3)$(4) :), 

+ (: $(l)$(2)$(3)n(4) :)„ (32) 

|(: n(l)$(2)$(3)$(4) :), = (: n(l)n(2)$(3)$(4) :), 

+ (: n(l)$(2)n(3)$(4) :), + (: n(l)$(2)$(3)n(4) :), 

+t(l)(: $(1)$(2)$(3)$(4) :), 

-3A[2((: $(1)$(2) :), + A(12))((: $(1)$(3) :), + A(13))((: $(1)$(4) :), + A(14)) 
+ (: $(1)$(1)$(2)$(3) :),((: $(1)$(4) :), + A(14)) 

+ (: $(1)$(1)$(2)$(4) :),((: $(1)$(3) :), + A(13)) 

+ (:$(1)$(1)$(2) $(1)$(3)$(4) :), 

+ (: $(1)$(1)$(3)$(4) :),((: $(1)$(2) :), + A(12)) 

+ (: $(1)$(1)$(3) $(1)$(2)$(4) :), + (: $(1)$(1)$(4) $(1)$(2)$(3) :), 

+ (: $(!)$(!) :),(:$(1)$(2)$(3)$(4) :),] 

-6A$o[((: $(1)$(2) :), + A(12))(: $(1)$(3)$(4) :), 

+ ((: $(1)$(3) :), + A(13))(: $(1)$(2)$(4) :), 

+ ((: $(1)$(4) :), + A(14))(: $(1)$(2)$(3) (33) 
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^(: n(l)n(2)$(3)$(4) :), = (: n(i)n(2)n(3)$(4):), + (: n(i)n(2)$(3)n(4):), 

+ (1 + V^ 2 )i{l){-. $(l)n(2)$(3)$(4) :), 

-3A(1 + V^2)[2{: $(l)n(2) :),((: $(1)$(3) :), + A(13))((: $(1)$(4) :), + A(14)) 

+ (: $(l)$(l)n(2)$(3) :),((: $(1)$(4) :), + A(14)) 

+ (: $(l)$(l)n(2)$(4) :),((: $(1)$(3) :), + A(13)) 

+ (: $(l)$(l)n(2) $(1)$(3)$(4) :), + (: $(1)$(1)$(3)$(4) $(l)n(2) :), 

+ (: $(1)$(1)$(3) $(l)n(2)$(4) :), + (: $(1)$(1)$(4) $(l)n(2)$(3) :), 

+ (:$( 1 )$( 1 ) $(l)n(2)$(3)$(4) :),] 

_6A$o(1 + V^2)[{-. $(l)n(2) $(1)$(3)$(4) :), 

+ ((: $( 1 )$( 3 ) :), + A(13))(:$(l)n(2)$(4) :), 

+ ((: $(1)$(4) :), + A(14))(: $(l)n(2)$(3) (34) 

|(: n(i)n(2)n(3)$(4):), = (: n(i)n(2)n(3)n(4):), 

+ (1 + Pi 2 + V^s)i{l){: $(l)n(2)n(3)$(4) :), 

-3A(1 + V ^2 + V^s)[2{: $(l)n(2) $(l)n(3) :),((: $(1)$(4) :), + A(14)) 

+ (: $(l)$(l)n(2)n(3) :),((: $(1)$(4) :), + A(14)) 

+ (: $(l)$(l)n(2)$(4) $(l)n(3) :), + (: $(l)$(l)n(2) $(l)n(3)$(4) :), 

+ (: $(l)$(l)n(3)$(4) $(l)n(2) :), + (: $(l)$(l)n(3) $(l)n(2)$(4) :), 

+ (:$(1)$(1)$(4) $(l)n(2)n(3) :), 

+ (:$( 1 )$( 1 ) $(l)n(2)n(3)$(4) :),] 

_6A$o(1 + V ^2 + V^s)[{: $(l)n(2) $(l)n(3)$(4) :), 

+ (: $(l)n(3) $(l)n(2)$(4) :), + ((: $(1)$(4) :), + A(14))(: $(l)n(2)n(3) :),] 

+ ^AA(1 - 3),5(2 - 3)[(: $(3)$(4) :), + A(34)], (35) 

4(: n(l)n(2)n(3)n(4) :), = (1 + V ^2 + V^s + Pi4)t(l)(: $(l)n(2)n(3)n(4) :), 

at 

-3A(1 + Vr 2 + ^3 + Vr,)[2{: $(l)n(2) $(l)n(3) $(l)n(4) :), 

+ (: $(l)$(l)n(2)n(3) $(l)n(4) :), + (: $(l)$(l)n(2)n(4) $(l)n(3) :), 

+ (: $(l)$(l)n(2) $(l)n(3)n(4) :), + (: $(l)$(l)n(3)n(4) $(l)n(2) :), 

+ (: $(l)$(l)n(3) $(l)n(2)n(4) :), + (: $(l)$(l)n(4) $(l)n(2)n(3) :), 

+ (: $(!)$(!) $(l)n(2)n(3)n(4) :),] 

-6A$o(l + ^2 + ^3 + $(l)n( 2 ) $(l)n(3)n(4) :), 

+ (: $(l)n(3) $(l)n(2)n(4) :), + (: $(l)n(4) $(l)n(2)n(3) :),] 
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CO I CM 


( 36 ) 


A[,5(l - 3).5(2 - 3)(: $(3)n(4) :), + ,5(1 - 4),5(3 - 4)(: n(2)$(4) :) 
+^(1 _ 4)^(2 - 4)(: n(3)$(4) :), + ,5(2 - 4),5(3 - 4)(: n(l)$(4) 
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